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Abstract. Dynamic models of M-type Mira variables predict the occurrence of water "shells", i.e. of zones of 
high H2O density and high H2O absorption inside the stellar atmosphere. The density, position and width of 
these shells is closely correlated with different types of two-component shapes of the intensity distribution on 
the disk in the H, K and L near-continuum bandpasses. We investigate these correlations and highlight the role 
of a spectrophotometric H2O index that warns against serious complications in diameter measurements in the 
case of substantial water contamination of the bandpass of observation. Simultaneous spectrophotometric and 
interferometric measurements may allow observers to estimate real continuum diameters more precisely. 
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1. Introduction 

Studies of Mira variables are of considerable astrophysi- 
cal interest as these objects are pulsating stars undergoing 
rapid mass loss and hence play an important role in the 
enrichment of the interstellar medium. Despite substan- 
tial progress in recent years, many fundamental aspects 
of the physics and structure of Miras still remain unclear. 
A summary of the present state of knowledge and of cru- 
cial problems that have to be solved was given, e.g., in 
a recent review of Scholz (2002 = S02). Among the most 
urgent problems is the determination of the stellar diam- 
eter that is needed for assigning to the star an effective 
temperature T e g, for placing the star in the HR diagram, 
and for understanding pulsation. 

Since the pioneer observations of monochromatic radii 
of o Cet and R Leo by Labeyrie and co-workers (Bonncau 
& Labeyrie 1973; Blazit et al. 1977; Labeyrie et al. 1977; 
Bonneau et al. 1982) and the first studies of dynamic 
Mira model atmospheres (Scholz & Takeda 1987; Bessell 
et al. 1989), it has been known that the atmospheres of 
M-type Mira variables are geometrically very extended 
configurations. Layers of formation of strong absorption 
features may be two or more times as distant from the 
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star's centre as the continuum-forming layers. For such an 
extended-atmosphere star, different diameter definitions 
may be given (Baschek et al. 1991) and the observed size 
of the stellar disk depends on wavelength. The most com- 
mon definition of "the" stellar radius, that describes the 
global dimensions of the star and is used in modelling the 
interior and atmospheric structure, is the Rosseland ra- 
dius R which marks the position of the layer at which the 
Rosseland optical depth equals unity, R — r(T-n oas =l). 

Clearly, the Rosseland radius is not an observable 
quantity, but R is usually close to the continuum optical- 
depth radius i? C ont — r(T cont =l) which depends only 
slightly on wavelength in the near-IR regime and which 
in turn may be related via model considerations to an ob- 
servable intensity radius (e.g. Bessell et al. 1996 = BSW96; 
Hofmann et al. 1998 = HSW98; Jacob & Scholz 2002 = 
JS02; S02). Substantial deviations of Rosseland radii from 
continuum radii are found in very cool M-type Miras in 
which the Rosseland extinction coefficient is noticeably af- 
fected by strong molecular absorption in high layers (e.g. 
HSW98; JS02; S02). The intensity radius (at continuum 
wavelengths or in an observer's filter passband) is de- 
fined in terms of the shape of the centre-to-limb variation 
(CLV). CLV reconstruction from interferometric data is 
in principle possible but not yet feasible with presently 
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available accuracies and baseline coverage. In practice, an 
intensity diameter is derived by fitting observed visibilities 
with model-predicted visibilities which are characterized 
by a radius- type parameter. Since molecular layers in the 
outer atmosphere significantly modify the pure-continuum 
CLVs, the corresponding visibilities and the resulting di- 
ameters, the choice of an appropriate atmospheric model 
is critical. 

Layers of high molecular density may often be consid- 
ered as a "molecular shell" which is more or less detached 
from continuum-forming layers depending on the temper- 
ature and density stratification of the Mira atmosphere. 
Indeed, attempts of interpreting near-IR CLVs in terms 
of H2O "shells" (e.g. Mennesson et al. 2002a, b) or visual 
light curves in terms of TiO "shells" (Reid & Goldston 
2002) are quite successful. Though this approach is useful 
and illustrative, one must keep in mind that these H2O or 
TiO shells are not circumstellar but are embedded in the 
stellar atmosphere, and quantitative interpretation of the 
effects upon CLVs and spectra may only be achieved in 
terms of molecule densities within the framework of the 
star's dynamic model atmosphere (e.g. Bessell et al. 1989; 
BSW96; HSW98; Hofner et al. 1998; Woitke et al. 1999). 
Model-predicted, shell-type molecular layers of water and 
various metal oxides were quoted by Scholz (1992; mod- 
els of Bessell et al. 1989) and Woitke et al. (1999). ISO 
and ground-based spectroscopic observations also point to 
warm layers of H2O (e.g. Hinkle & Barnes 1979; Tsuji et 
al. 1997; Yamamura et al. 1999; Matsuura et al. 2002) and 
other molecules. 

In this paper we focus on the characteristics of 
CLV shapes in the water-contaminated continuum band- 
passes in the near-IR. We concentrate our discussion on 
the K bandpass (HSW98; rectangular filter with A c = 
2.195/ito, AA = 0.40/xm) as this is the filter which is most 
commonly used for interferometric observations of Mira 
variables. We also relate the different classes of CLVs pre- 
dicted by the models to the time-dependent structure of a 
pulsating star and to observable spectrophotometric prop- 
erties. A simple classification is described in Sect. [21 In 
Sectional we demonstrate that shells of high water vapour 
densities are found in the outer layers of most models as a 
natural consequence of the density and temperature strat- 
ification. Their location and extent is shown to relate to 
the shape of the CLVs, establishing water vapour as an 
important actor on the scene of angular diameter mea- 
surements. In Sect. 0J we show that spectrophotometric 
indices sensitive to water bands can in principle warn the 
interferometric observer against the presence of visibility 
distortions and provide estimated corrections to angular 
diameter values derived in the uniform disk approxima- 
tion. However, the reliability of these corrections depends 
on the precision and the parameter coverage of predicted 
model spectra for individual observed stars. In a forthcom- 
ing paper (Tej et al. 2002 in preparation) , we will provide 
an extensive comparison with empirical optical to near-IR 
spectra. 



2. Classification of CLVs 

In this paper we investigate the CLVs predicted by the 
models of BSW96 and HSW98. For the HSW98 models 
(P, M, O series), the pulsation is driven by complete self- 
excited configurations whereas, for the BSW96 models (Z, 
D, E series), a conventional piston approach to the sub- 
atmospheric layers is used. The density stratifications re- 
sults from the shock front driven outflow and the subse- 
quent infall of matter. Non-grey temperature stratifica- 
tions are computed in the approximation of local ther- 
modynamic and radiative equilibrium. Computationally 
smeared-out density discontinuities occur at shock front 
positions. We refer the reader to BSW96 and HSW98 for 
detailed discussions on the construction of these dynami- 
cal models. 

The models of Bessell et al. (1989), BSW96 and 
HSW98 predict a large variety of CLV shapes. These 
range from nearly uniform and moderately darkened disks 
to complex CLVs with Gaussian-like or two-component 
appearance. The model study of JS02 shows that pure- 
continuum intensity diameters in the near-IR may be ac- 
curately defined in terms of the position of the steep flank 
of the CLV. The resulting visibilities of interferometric 
observations of the stellar disk may readily be fitted by a 
moderately darkened disk, and simple uniform-disk (UD) 
fits lead to reliable continuum diameters with errors of the 
order of only a few percent. This study also shows, how- 
ever, that the standard near-IR continuum bandpasses in 
the 1.2 to 4.0/im range (J, H, K, L) are noticeably con- 
taminated by atmospheric water absorption in Miras with 
low effective temperatures and strongly extended atmo- 
spheres. The CLV has a two-component shape consist- 
ing of a moderately darkened inner continuum portion 
and an outer tail- or protrusion-type extension. The two- 
component shape of the brightness profile generated in 
the Rayleigh- Jeans part of the Planck function by mod- 
est molecular absorption can well be understood in terms 
of basic radiation transport (Scholz 2001). The result- 
ing visibility appears distorted as compared to the pure- 
continuum visibility. Pronounced protrusions may mimic 
a Gaussian-like shape of the full CLV. The detailed dis- 
cussion of JS02 shows that these effects lead to serious 
consequences for measuring continuum diameters. Here, 
we restrict our investigations and discussions to these two- 
component CLVs. Figure H displays the model CLVs for 
the K passband. The parameter of the non-pulsating "par- 
ent" star and the time series of the Mira models are dis- 
cussed in HSW98 and quantities relevant to our study are 
listed in TablcTT] and TableEJ The CLVs can be broadly 
classified into three groups based on visual inspection. 
Note that this is a purely qualitative description ; a more 
quantitative one will be introduced in Sect.0] 

— Tail Type (T): The intensity of this tail is model depen- 
dent. The turn-over from the flank to the tail feature 
is usually located near or below 0.2 in the normalized 
intensity scale. The angle of turn-over also depends on 
the model. For instance for P20 the turn-over is sharp 
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Table 2. Parameters of the time series of the Mira models (adapted from BSW96 & HSW98) and classification of 
model CLVs. The columns: visual phase <fi v i s ; luminosity L; Rosseland radius R; 1.04 near-continuum radius i?i.o4 
(JS02); effective temperatures T e g(R) and T e g (Rim)] classification and remarks. An abbreviated nomenclature (letter 
defining the series followed by the pulsation phase (x 10)) is adopted for the model names. 
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Note: P(?) indicates protrusions with atypical characteristics (see text). N(?) indicates quasi-normal CLVs, that do however 
display a small pedestal. 



Table 1. Properties of the Mira "parent" star series 
(taken from BSW96 & HSW98). The columns: pulsation 
mode - fundamental (f) or overtone (o); period P; parent 
star mass M; luminosity L; Rosseland radius R p ; effective 
temperature T ff oc (L/i? p 2 ) 1 / 4 . 
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but for E10 it is comparatively smoother. The extent 
of the tail can also vary as is seen for P40, where the 
tail is shorter compared to P20. 

— Protrusion Type (P): When the strength of the tail 
increases it appears as a protrusion of the inner disk. 
The protrusion appears at a high intensity level. The 
nature of this feature is also strongly model dependent. 
It can either be in the form of a 'bump' as in P35 or a 
smoother feature as in M15. 

— Normal type (N): These display neither tails nor pro- 
trusions and are not two-component CLVs. They could 
be modelled by moderately darkened disks. 

Table 2 lists the class assigned to each K band CLV, 
following the above definitions. The strength of the sec- 
ond CLV component indicates the amount of water con- 
tamination and is strongly wavelength-dependent, but the 
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Fig. 1. The model CLVs. The radius is given in arbitrary units (consistent within each panel). 



basic characteristics of the brightness profile are the same 
in the H and L bandpasses with less and more water, re- 
spectively, than the K band (JS02). In the near-maximum 
models of the M series, a protrusion-type attachment to 
the inner CLV is seen which is inconspicuous and almost 
merging with the inner CLV, in particular in the H band 
where water absorption is quite modest. In the very cool 
first-overtone O models, the H 2 O-generated protrusion is 
strong, but it appears only in the near-minimum model 
O08 as a clearly separable component wheareas it blends 
with the inner CLV in O05 and O10. 



It is evident from the plots and the table that a single 
parameter, such as the phase, alone does not explain the 
presence of tails or protrusions. For instance, the near- 
minimum models of the relatively warm Z and D series 
have tails whereas the cooler P models have protrusions 
near-minimum light. Average temperature and pulsation 
mode also play a significant role. If, for the purpose of 
ordering the classes, we call protrusions stronger features 
than tails, there is a clear correlation between the strength 
of the feature and the phase within individual model se- 
ries. The strongest features are found near-minimum. Note 
that for the overtone series, which have larger radii and 
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lower T c g than fundamental pulsators, models show a 
global shift of the minimum luminosity to ~ 0.8 (BSW96; 
HSW98). In the subsequent sections it will become evi- 
dent that the diversity of the CLV shapes seen is due to a 
combination of the shock propagation, the shape and ex- 
tent of the water shells and the density and temperature 
stratification. 

3. Water "shells" and their implication on CLV 
shapes 

Recent interferometric observations of Miras by 
Mennesson et al. (2002b) show a strong increase in 
the measured UD diameter between the K' and the L' 
bands. These authors propose a two-layer scenario with a 
classical photosphere and a detached and extended (~ 3 
stellar radii) gas layer to account for this chromatic varia- 
tion in size. Similar two-layer models are also invoked by 
Yamamura et al. (1999) to model the ISO/SWS spectra 
of o Cet and Z Cas. A two-layer "slab" model with plane 
parallel uniform molecular layers of H2O placed one 
above the other is shown to accurately fit the observed 
spectrum. 

We investigate the presence of these water shells in 
time series of Mira models and their implication on the 
shape of the CLVs. To display the shells, we use the quan- 
tity p x Pwater / Pgas, where P gas is the total gas pressure 
and P W ater the partial pressure of H2O. It measures the 
density of water molecules. In Figures |21 and |21 we plot 
the P models illustrating the protrusions and tails respec- 
tively. For Figure |2 we show the P05, P15 and the P35 
model CLVs (P25 is essentially identical to P05) and for 
FigureEl we plot the P10, P20 and P40 models (P30 is 
similar to P40 with a slightly broader continuum disk). 
The upper panels show this shell parameter as a func- 
tion of the radius, the middle panels plot the respective 
CLVs and the lower panels show the density structure for 
these models. Note that the computational smoothing of 
the shock fronts slightly affects the detailed behaviour of 
the corresponding areas in all panels. The profile of the 
water shells clearly plays a role in defining the shape and 
strengths of the tail/protrusion component of the CLVs. 

For the protrusion-type P models, which occur dur- 
ing minima, the effective temperature is low and there 
is a steep density gradient after the position of the in- 
nermost shock front. Shells are formed relatively close to 
the continuum-forming layers. For the P05 model, there 
is a narrow "peak-shell" and the density of the water 
molecules goes down abruptly at the density discontinuity 
position of the outer shock front. This manifests itself as a 
sharp cut-off point on the protrusion feature beyond which 
the intensity is small. Similar behavior is seen for the P35 
model where the density of water molecules sharply de- 
creases at the position of the outer shock front. The wa- 
ter shell in this case is much broader and of low density 
as compared to the P05 model. This is reflected in the 
CLV shape where the water component of the CLV de- 
creases rather slowly. In the P15 model, there is no outer 



shock front, hence, the water shell is much more extended. 
The above discussion suggests that the extension of these 
shells depends significantly on the nature of propagation 
of shocks. 

The water shells for the tail- type P models are of differ- 
ent nature. These occur during maxima and for higher ef- 
fective temperature. Tails correspond to shells with signif- 
icantly lower water density than protrusion. The shells are 
detached from the continuum-forming layers and are lo- 
cated in higher atmospheric regions. P10 has an extended 
but very-low-density shell. P20 has a higher density and 
there is an abrupt cut-off in the shell coinciding with the 
position of the outer shock front. The corresponding CLV 
has a constant tail feature up to the cut-off point. The shell 
is much more extended for the P40 model and without the 
occurrence of a outer shock, the CLV shape remains con- 
stant out to larger radii before fading to zero. 

In Figures 01 and [3 we show similar plots for the high- 
mass overtone O and fundamental M models (M15 is simi- 
lar to M05 and hence not plotted). As can be seen from all 
these plots there are a variety of shapes, strengths and ex- 
tensions of these water shells which agree quite well with 
the variety of shapes seen in the CLVs. A point worth men- 
tioning here is the presence of a sharp spike seen in the 
shells (for e.g. in the P40, O08, M05 models) close to the 
continuum-forming layers. This is due to the enormous in- 
ward increase of p while P wa ter / Pgas is still non-zero. This 
does not affect the two-component CLVs but results in 
some additional darkening of the normal continuum disk. 

Apart from the very good one-to-one correlation seen 
between the water shells and the shapes of the CLVs, 
in general we see that the shells are detached and more 
extended for the tail-type models as compared to the 
protrusion-type ones. This is in gross agreement with the 
results of Yamamura et al. (1999) who surmise that H2O 
layers are more extended at maxima. Mennesson et al. 
(2002b) derive a photospheric radius of 10 mas and a H2O 
envelope extending from 15 mas to 27 mas for R Leo. This 
"gap" of 1.5 Ri between the continuum layer and the H2O 
envelope is in good agreement with the shells seen for the 
tail-type P models. Occurrence of such "gaps" is the nat- 
ural consequence of the temperature-density requirement 
of the equation of state for H2O formation. 

The density stratification of various molecules exhibits 
shell-like structures at locations similar to those of H2O, 
though with differing relative strengths and extents (cf. 
Woitke et al. 1999). For the models used here, JS02 found 
that molecular contamination of the H, K and L* con- 
tinuum bandpasses is essentially due to H2O with small 
contributions of OH and CO in the H bandpass. Hence, 
molecules other than H2O probably do not appreciably 
affect the CLVs and the corresponding visibilities in the 
near-IR passbands. However, one should be aware that 
different filter profiles may result in different molecular 
contributions, that the parameter range and phase-cycle 
coverage of the models is limited, and that model assump- 
tions and adopted model opacities may play a role. For 
instance, the BSW96 and HSW98 models used here as- 
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Fig. 2. The relation between the the water "shells" (up- 
per panel), the CLV shapes (middle panel) and the den- 
sity structure (lower panel) for the near-minimum P model 
CLVs. The line identifications are P05 - solid; P15 - dot- 
ted; P35 - dash-dotted. R p is the radius of the static par- 
ent star (Table 1). In the middle panel the positions of the 
continuum-forming layer (-R1.04) are shown with crosses. 
In the lower panel the location of the shock fronts (com- 
putationally smeared out) are marked with filled circles. 



sume local thermodynamic equilibrium, whereas Woitkc 
et al. (1999) have shown significant deviations from equi- 
librium for four molecules including H2O. Also, variations 
of metallicities and of C/O ratios (transition to MS- and 
S-type Miras) are not considered in the present models. 

Finally, some interferometric and spectroscopic obser- 
vations indicate that inner edges of dust shells might 
be as close as 2 or 3 continuum radii from the star's 
centre (e.g. Danchi et al. 1994; Danchi & Bester 1995; 
Lobel et al. 2000; Lorenz-Martins & Pompeia 2000). M- 
type Mira models with thermodynamically coupled dust 
(Bedding et al. 2001) show two-component CLVs simi- 
lar to those produced by molecular layers for the 1.04 /im 
bandpass, a bandpass that is not significantly contami- 
nated by molecules otherwise. The distortions are shown 
to be pronounced around and below 1 /im but the effects 
are almost negligible for the H, K and L* passbands. The 
situation is likely to depend strongly on the treatment of 
dust condensation, and further theoretical work will be 
useful. For C-type Mira atmospheres, dust condensation 
is significant in existing models (e.g. Hofncr et al. 1998), 
but these stars lie outside the scope of our paper. 
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Fig. 3. Same as in Fig. [21 for the near-maximum models. 
P10 (solid), P20 (dotted) and P40 (dash-dotted). 
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Fig. 4. Same as in Fig.|21 for high-mass overtone models 
O05 (solid); O08 (dotted) and O10 (dash-dotted). 

4. Spectral indications of visibility distortions 

In the previous section we have discussed the occurrence 
of water shells, due to absorption in the outer atmospheric 
layers, and their effects on the shape of the CLVs. JS02 
show that this molecular-band contamination can dis- 
tort visibility curves and make the measurement of pure- 
continuum diameters of Miras very difficult. Visibility dis- 
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Fig. 5. Same as in Fig.[21for high-mass fundamental mod- 
els M05 (solid); M10 (dotted) and M20 (dash-dotted). 



tortions are not seen for normal CLVs or those with tiny 
tail-like features (e.g. P10). In such cases, the diameter 
estimated by fitting a UD or fully darkened disk to the 
measured visibilities would represent the near-continuum 
value. But for models which have tails of high intensity 
(e.g. P20) the visibilities are highly distorted. The dis- 
tortions are even more pronounced and are more or less 
Gaussian in shape for models with protrusion-type CLVs 
(e.g. P35). Consequently, the fit diameter (by means of 
a UD or a fully darkened disk) is appreciably depen- 
dent on the spatial frequency (i.e. the baseline of obser- 
vation). Observationally, this effect is seen in the shape 
of high-precision visibilities with suitable baseline cov- 
erage (e.g. Perrin et al. 1999; Weigelt et al. 2002) or 
in the baseline dependence of UD fit diameters, and in 
the bandpass-dependence of near-continuum UD fit diam- 
eters (e.g. Tuthill et al. 2000; Mennesson et al. 2002b; 
Thompson et al. 2002; conference reports listed in JS02). 

Clearly, the choice of an appropriate atmospheric 
model and of the corresponding CLV is a critical step 
in the determination of a stellar diameter. One way of 
constraining this choice is to combine interferometric ob- 
servations with many baselines or many wavelengths, and 
to search a set of atmosphere models for the best fitting 
ones. With today's instrumentation and relatively poor 
parameter coverage of the models, this remains cumber- 
some or impossible. Combining a few interferometric mea- 
surements with simultaneous spectroscopy should turn out 
to be easier. This justifies the search for spectrophotomet- 
ric indicators of distorted visibilities. 

As it now appears, water molecules are the main cul- 
prit of the two-component nature of the near-IR CLVs. 
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Fig. 6. Strength of the tail/protrusion feature as a func- 
tion of the H2O index for the H, K and L* bandpasses. 
The symbols are: filled triangles - normal type CLVs; open 
circles - tail-type CLVs and filled circles - protrusion-type 
CLVs. The 4 models with atypical P(?) protrusions in 
Table 2 (see text) are marked by arrows: M10 (H2O index 
0.84), M20 (1.06), O05 (1.83), O10 (1.45). 



Hence, we look for spectrophotometric features which 
would quantify the amount of water present in Miras. The 
H2O index defined by Persson et al. (1977) and Aaronson 
et al. (1978) is a good indicator, that has been frequently 
used for classification purposes (Elias et al. 1982; Terndrup 
et al. 1991; Langon & Rocca-Volmerange 1992). It mea- 
sures the ratio of the fluxes measured through narrow- 
band filters centered, respectively, in the wing of the 
1.9 /urn water band (A = 2.0 /im) and in the K band 
(A = 2.2 /im). The index is expressed in magnitudes 
and is zero for Vega. To quantify the strength of the 
tail/protrusion feature of the CLVs, we divide the CLV 
into two parts - the inner disk, and the tail- or protrusion- 
type extension. The ratio of the area of the extension with 
respect to the inner disk gives an estimate of the feature 
strength. 

In Figure El we show the feature strength as a func- 
tion of the H2O index for the H band (HSW98; rectan- 
gular filter, A c = 1.633^m, AA = 0.30^m) in the upper 
panel, the K band in the middle panel and the L* band 
(HSW98; rectangular filter, A c = 3.799/im, AA = 0.60/^m) 
in the lower panel. This figure shows the direct correlation 
between the strength of the CLV extensions and the H2O 
index. The strength increases with the increase in the H2O 
index. For the H and K band CLVs the ratio of the areas 
increases slowly from the normal-types to the tail-types 
and then there is a steep rise for the protrusion-type exten- 
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sions. The three large-ordinate points marked by arrows 
in the upper panel of Figure 6 are atypical P(?) protru- 
sions in Table 2, which are particularly hard to separate in 
an objective manner from the inner CLV for the H band 
(see Section 2). Water contamination in the near-IR band- 
passes is stronger in the L* band as compared to the H 
and K bands (JS02). This manifests in increased strength 
of the tail/protrusion feature of the CLVs, as clearly seen 
in the lower panel of Figure 

The spectrophotometry index of H2O gives a satisfac- 
tory estimate of the strength of the tail/protrusion feature 
of the CLV. However, the CLVs are not directly observ- 
able quantities but the shapes have a direct implication 
on the measured visibilities. JS02 have done a detailed 
study of the effects of phase, cycle and baseline on the 
measured continuum diameters of different Mira models. 
They fit UD visibility to the full model-predicted visibili- 
ties by means of a least square routine and also carry out 
single point fits of the UD to the model visibility at three 
spatial frequencies the choice of which covers the entire 
range of various model visibilities. The least square fit- 
ting to the full model takes care of the arbitrary choice 
of baselines and the three point fitting gives the baseline 
dependence of fitted diameters. The near-IR bands show 
substantial deviation of the fitted UD diameters from the 
near-continuum values. Apart from uncertainties of model 
treatment of molecular and dust opacities at very low 
temperatures, i?i.o4 (the monochromatic radius defined 
at T1.04 = 1 in the 1.04 /im continuum window) is very 
close to the continuum value. We refer the reader to JS02 
for a detailed discussion of the points outlined above. As 
an obvious next step, we explore the correlation between 
the H 2 index and the deviations of fitted UD diameters 
from the continuum values. We use the K band results of 
Figures 8, 9 & 10 of JS02 for our study. 

In the lower panel of Figure we plot the ratio of the 
least squares fit UD radius (-Risq) to the near-continuum 
radius (-R1.04) as a function of the H2O index. Although 
the dispersion is large, the deviations of the fitted UD 
diameters from the continuum values distinctly increases 
with water contamination. The relation tightens if the 
overtone models (O and E series) are left out. A linear 
least squares fit to the model points excluding the over- 
tone models yields 

iW-Ri.M = 0.19 H 2 + 0.89 (1) 

In the upper panel we plot the correction factor for long 
baseline observations {R\ on g/ Ri.qa)- This baseline is po- 
sitioned near the first null of the visibility function. It is 
more sensitive to the details of the CLV structure. The 
least squares fit to this gives 

R\onjR\.QA = 0.11 H 2 + 0.90 (2) 

Similar trends are seen for the medium and short baselines 
(not plotted). The ratios of interferometric UD radii to the 
Rosseland radius R (instead of Rim) show much more 
scattered relations to the H 2 index. This is expected be- 
cause of the nature of the Rosseland opacity which is a 
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H 8 Index 

Fig. 7. In this figure we plot the deviations of the UD fit- 
ted diameters from the the -R1.04 for the K band. The data 
points are taken from Figures 8,9 & 10 of JS02. The lower 
panel shows the results from the least squares and the 
upper panel shows the long baseline values. The symbols 
are P models - filled circles; D models - filled triangles; M 
models - filled squares; O models - open circles; E models 
- open triangles 

harmonic-type mean opacity and is sensitive to extended 
regions of strong molecular extinction which tends to close 
the near-continuum windows, to increase the Rosseland 
opacity and, hence, to increase the Rosseland radius as 
compared to the continuum radius. 

The above two linear fits emphasize the fact that a 
simultaneous spectral observation of water contamination 
could yield a correction factor for deriving continuum di- 
ameters from the UD fitted values. A couple of caveats 
should be kept in mind. First, the dispersions in the plots 
remain large. The situation should improve when more 
models become available, allowing us to explore the origin 
of the spreads. For instance, the correlation may depend 
on other measurable quantities, such as pulsation prop- 
erties, as suggested by the deviant behaviour of the O 
series. Second, the exact shape and depth of water bands 
in spectra is difficult to compute theoretically, in partic- 
ular in extended, dynamical atmospheres. The reliability 
of the numerical correction factors given above depends 
on the accuracy the model spectra and on the range of 
model parameters explored. As will be shown in a forth- 
coming paper, current models are successful in reproduc- 
ing the global features of empirical spectra of Miras and 
the range of observed water indices, but uncertainties re- 
main large in the time-dependent depths and shapes of in- 
dividual bands. Once fully satisfactory models exist, other 
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spectrophotometric measures of the water bands may turn 
out to be more sensitive to the radius ratios of Fig.0 The 
series of BSW96 and HSW98 display a tendency for a 
broader 1.9 /im band in models with protrusions than in 
models with tails that remains to be studied in detail. 

5. Conclusion 

In this paper we have classified the two-component near- 
IR CLVs predicted by the pulsating Mira models of 
BSW96 and HSW98 based on a qualitative visual inspec- 
tion. The models are shown to naturally produce water- 
rich shells depending on temperature and density strat- 
ification of the atmospheric models. The propagation of 
shock fronts plays a decisive role in the structure and ex- 
tent of the shells and has direct effect on the shapes of the 
CLVs. The model-predicted strength of the H2O index, 
quantifying the 1.9 /im water feature, gives a good indica- 
tion of the nature of the CLV in terms of the strength of 
the tail/protrusion features with respect to the inner con- 
tinuum disk. The increase in water contamination from 
the K band to the L band is evident from the increase in 
the feature strength for the L band CLVs. 

Using the results of JS02, we show that the deviations 
of the fitted UD diameters from the continuum values are 
correlated with the H2O index, though with a significant 
dispersion. The results are based on a moderate number 
of models covering a limited range of stellar parameters, 
phases and cycles, which prevents us from identifying the 
origin of the scatter. The exact form of this correlation still 
has to be determined from future high-precision visibility 
observations, and the usefulness of a simple linear approx- 
imation as given in Eqs. (1) and (2) still has to be checked 
both empirically and with more models. Until then, simul- 
taneous spectroscopic observations of water bands may 
at least warn the interferometric observer and provide a 
rough estimate of corrections to be applied to the UD fit 
diameters. 
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